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I. INTRODUCTION 

A. Preamble 

The network study is an old field of research, which 
has recently become particularly active (see [1]). This 
is largely due to the opening of an access to rich data 
on natural systems: the World-Wide Web, the Internet, 
the various biological networks (gene transcription, cell 
metabolism, protein interactions), the sociological net- 
works (citation index, collaborations, phone calls), etc. 
Interesting empirical regularities have been observed, like 
the small-world property of networks or the frequently 
observed scale-free nature of degree distributions. On 
the theory side, the natural conceptual framework for 
network research is the graph theory: the nodes and the 
links of a graph represent the active agents and the inter- 
actions, respectively. However, it was soon realized that 
the classical graph theory [2] is inadequate and has to 
be generalized, in order to cover the new reality. This 
observation has triggered, in turn, an intense theoretical 
activity, that led to the construction of a number of in- 
sightful models. It is clear, however, that much remains 
to be done. 

As recalled in ref. [3], in studying complex systems 
one can adopt one of the two complementary approaches: 
the diachronic and the synchronic one. In the former, 
one focuses on the time evolution of the system, which 
helps discovering the dynamics at work. In the latter, one 
considers an ensemble of dynamically similar systems at 
a fixed large time, which helps identifying the generic 
structural traits. 

Most of the recently constructed models adopt the di- 
achronic approach and deal with growing networks [4]- 
[20]. This led, among others, to the remarkable discovery 



of the preferential attachment rule and of its crucial role 
in network evolution [4, 5]. The mathematics of the di- 
achronic approach rests largely on master equations and 
related manipulations. The synchronic approach has also 
been advocated and various static statistical ensembles 
have been constructed and studied [3],[21]-[30]. Here, 
the mathematics is that of the equilibrium statistical me- 
chanics and probability theory. One of the aims of this 
paper is to help establishing a bridge between the two sets 
of results. It will be seen that the widely accepted dis- 
tinction between growing and equilibrium networks is not 
really correct. The truly distinctive property of growing 
networks, their causal structure, can be usefully incorpo- 
rated in a static model. We limit ourselves to tree graphs, 
because in this case the use of the methods of analytic 
combinatorics enables one to get exact solutions. 

We shall discuss at length the geometry of connected 
tree graphs endowed with a causal structure. Our discus- 
sion will cover, in particular, some of the recently popular 
growing network models [4, 6] , but we shall proceed con- 
sidering static statistical ensembles of trees, employing a 
formalism which proved to be useful in deriving generic 
results for other random geometries. For example, we 
shall show that the causal trees have the "small world" 
property: their Hausdorff dimension is infinite. The ad- 
vantage of our approach is that is enables one to show at 
once that this property holds in a wide class of models. 

The plan of the paper is as follows: In the next two 
subsections we recall some basic concepts and we define 
the observables to be calculated later on. In sect. II 
we begin by deriving a number of general results that 
hold for causal trees: The central role is played by 
a recursion relation satisfied by the partition function 
(sect. HA). This recursion relation enables one to find 
closed expressions for the degree distribution (sect. IIB), 
the ancestor-descendant correlation (sect. IIC) and the 
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two-point function measuring the distribution of shortest 
path lengths between pairs of nodes (sect. IID). In sect. 
HE we establish contact between our formalism and the 
growing network models. This is used in sect. IIF, where 
we calculate explicitly a set of observables for some sim- 
ple growing networks. In sect. Ill we compare the causal 
model with the maximally random one: In sect. IIIA we 
collect all general formulae. In sect. IIIB we calculate 
the observables, assuming the same microstate weights 
as in sect. IIF, and we find dramatically different results. 
This illustrates the importance of graph symmetries. In 
sect. IIIC we show how to construct a maximally ran- 
dom model with the same degree distribution as a given 
causal model. We stress, however, that the identity of 
degree distributions does not imply that geometries are 
similar. In particular, the average shortest path lengths 
are very different: generically they scale like a power of 
the number of nodes N in maximally random and like 
log N in causal models, respectively. We conclude briefly 
in sect. IV. 



B. Basic concepts 

Let us recall some definitions: A rooted tree is a tree 
with one marked node. A planted tree is a rooted one 
with an extra link attached to the root, so that the degree 
of the root is increased by unity. The other end of this 
extra link is not counted as a node, in a sense this end 
remains "free". The different tree ensembles are simply 
related and choosing to work with one of them is a matter 
of convenience. In this paper we deal with planted tree 
graphs [31]. 

By assumption, a label is attached to each node and 
two graphs with identical topology but labeled differently 
are considered different. We say a tree is endowed with a 
causal structure when the labels always appear in growing 
numerical order as one moves along the tree from the root 
towards an arbitrary node. These are the tree graphs we 
are interested in. 

Introduce a statistical ensemble of these trees. Let us 
denote by T a given topology and by L(T) the number of 
distinct causal labelings of T. We attach the same weight 

p(T) = w(m,n 2 ,...,n N ) (1) 

to each acceptable labeling. Here N is the total number 
of nodes, denotes the degree of the node i and w is 
some appropriate positive-definite function. A model is 
defined by choosing a particular form of w. 

It will be seen that the presence of a causal structure 
generates non-trivial observable inter-node correlations. 
Hence, it is of interest to discuss models where these spe- 
cific correlations do not interfere with correlations of a 
different origin. With this motivation, we assume in this 
paper that p(T) factorizes [32] 

N 

p(T) = Y[q ni (2) 



On the other hand, we keep q n as general as possible. 

The partition functions of canonical and grand- 
canonical ensembles arc defined in the usual way by sum- 
ming the weights of all possible microstates. Thus by 
definition : 

^ = ^X L{T)p{T) (3) 

T 

and 

Z{p)=Y,ZNe-» N (4) 

N 

The prefactor 1/N\ in cq. (3) is compensated by the 
number of terms in the summand. Indeed, the number 
of labeled causal trees with N nodes is (N — 1)! [33]. 
It is not difficult to convince oneself that in all cases of 
physical interest Zn grows exponentially with N, up to a 
power prefactor : 

z N ~ e i* N for N -» oo (5) 

and Z{p) develops a singularity at p = p (see sect. IIA 
for a more rigorous argument). One is primarily inter- 
ested in the regime controlled by the singular part of 
Z(p). Indeed, as Sp = p — p tends to zero one becomes 
increasingly sensitive to the behavior of trees with arbi- 
trarily large N. 

Of course, eq. (4) can be inverted by Laplace transfor- 
mation 

p,N rioo 

Z ^ = ^T- d <^ Z(Sfi) e s » N (6) 
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with the integration contour passing on the right of the 
singularity at dp — 0. It will often be convenient to work 
in the grand-canonical ensemble and to Laplace trans- 
form the result to the physically more interesting canon- 
ical ensemble at the very end. In practice, we shall always 
assume that N — > oo and we shall keep the leading term 
only. 

Notice, that there is an analogy between the labelings 
of graphs and the positions of a system in a (discrete) 
phase-space. Hence, p(T) and L(T) are the analogues 
of the weight of a microstate and of the corresponding 
phase-space volume, respectively. 

Little is specific to causal trees in the content of this 
section. One could repeat verbatim the above definitions 
and keep the same microstate weights in the context of 
a different tree ensemble. However, in this new ensem- 
ble the trees would have different symmetries and, there- 
fore, L(T) would be in general different. Consequently, 
the physics would be different too. We shall see later, 
comparing causal and maximally random trees, that in- 
troducing a new symmetry can change dramatically the 
geometry of generic graphs. 



C. Observables 

Let us define the observables we shall calculate in this 
paper : 

The most popular observable is the degree distribution. 
When the weight of a microstate has the factorized form, 
as in (2) , the degree distribution ir n is given by the simple 
and obvious formula 



N-\ n 



<91n zm 
dq n 



(7) 



The factor TV -1 above is included to have the distribution 
normalized to unity. 

In the thermodynamic limit N — > oo 



d/j, 
dq„ 



(8) 



The next observable is the correlation between node 
degrees, say the probability that a node has degree k 
when its neighbor's degree is I. In a causal tree, one of 
these nodes is an ancestor and the other a descendant. 

A very interesting observable is the Hausdorff dimen- 
sion dn controlling the scaling with N of the linear size 
of a typical tree : 



(9) 



One usually takes for r the distance between an arbi- 
trary pair of nodes. In rooted trees it is more natural to 
consider the distance separating a randomly chosen node 
from the root. One first calculates a specific two-point 
function C (r, fi) , the grand-canonical weight of all trees 
with a node separated from the root by r steps. Using 
C(r, /«) one finds 



Jo°° dr CV, M) 



(10) 



This quantity usually diverges when 5(j, — > 0. The behav- 
ior on the r.h.s. (right-hand side) of (9) is determined by 
observing that Sfi scales like TV -1 (see (6)). 



II. GEOMETRY OF TREES ENDOWED WITH 
A CAUSAL STRUCTURE 

A. Recursion relation 

We start by deriving a recursion relation for the parti- 
tion function zn- To this end we construct a new planted 
tree by attaching the k "free link-ends" of the planted 
trees Ti, . . . ,T k to a new root. We denote the resulting 
compound tree by T — Ti © • • • © T k (this is illustrated 
in FIG. 1 for k = 3). 

The number of ordered labelings of the compound tree 



is 



L(Ti © • • • © T k ) 



N\ 



1 



mi-'-Nklkl 



L(T 1 )---L(T k ) (11) 





FIG. 1: Construction of the compound tree T 1 ®T 2 @T 3 . The 
three old roots are connected to the new root with degree 
four. Notice that all the graphs are planted trees. The trees 
are not planar and therefore the relative position of branches 
in the compound tree is irrelevant. 



Here Ni denotes the total number of nodes in the tree 
Ti and N = Y,i N i- 0nc nas to g ive N + 1 labels to 
the nodes of the compound tree. However, the smallest 
label must be attached to the root. The remaining TV 
labels are arbitrarily distributed among the trees. This 
is the origin of the the multinomial factor. Permuting the 
trees Ti does not change the compound tree [35]. This 
explains the presence of the factor 1/kl. It will be seen 
that the composition rule (11) is a very basic property of 
the model. 

Obviously the weight of the new tree factorizes 

p(Ti © • • • © T fc ) = q k +ip{Ti) ■ ■ ■ p(T k ) (12) 

The partition function zn+i can be constructed by 
summing the trees of size smaller or equal to N : 



1 



ZN+l 



{N + iy. 



N! + ---+N k ,N 



fe=lTi,...,T fc 

x L(Ti © • • • © T fc )p(Ti © • • • © T fe ) 



(13) 



Inserting (11) into (13) and rearranging terms in the 
sum we obtain after trivial algebra : 

1 ..^ k 

Z N+1 = iv + 1 E E S N! + -+N k ,N Hz Ni (14) 



Adding z\ — q\ and summing both sides of equation (14) 
we get 



N *Ne 



Nn 



or finally 



where 



Z\n) = -e-»F{Z) 



qk 



Equation (16) can be integrated to give 



Jo 



da; 
F(x) 



(15) 



(16) 



(17) 



(18) 
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The function G(Z) is a positive monotonically growing 
function of Z, bounded from above (one can ignore the 
trivial case where all q n except qi and q 2 are zero). Hence 
[i is bounded from below: Z(fi) has a singularity at some 
yu = p,. Denote by x the radius of convergence of the 
series F(Z). The critical value of ji is given by 



M = -logG(x) 



(19) 



This formula holds also when the radius of convergence 
x is infinite, since all terms in the series (17) are posi- 
tive and the integral in (19) is convergent in all cases of 
interest : G(oo) < oo. 



B. Degree distribution 

The degree distribution calculated using (8) is 
1 q n f x dx 



7T„ 



G(x) (n-l)\J F(x) 



(20) 



Again, this formula is also valid when x = oo. 

Summing over n and using the definitions of _F and 
G one easily checks that 7r„ is normalized to unity, as it 
should. One further finds 



E 



n tt„ 



F{x)G{x) 



(21) 



On a tree, the r.h.s. should equal 2 . This is the case 
when F{x) diverges at x = x. Otherwise one encoun- 
ters a pathology, which looks similar to that appearing 
in some maximally random tree models (and in the so- 
called balls-in-boxes model, see [3, 34]), where working in 
the large N limit one misses singular nodc(s) contribut- 
ing term(s) of the type 7V _1 (5(n — cN). In the limit such 
nonuniformly behaving terms do not contribute to the 
normalization but do contribute to the r.h.s. of (21). It 
will be shown later that when F(x) < oo the average 
distance between nodes is finite. This means that singu- 
lar node(s) - with unbounded connectivity - are indeed 
expected to show up. In deriving (20) the large iV limit 
has been implicitly used and it is an educated guess that 
one again misses the singular node(s). 



C. Ancestor-descendant correlation 

Now, wc turn to the calculation of the ancestor- 
descendant degree correlation. It is obvious that an an- 
cestor plays the role of the root of the subgraph involving 
all its descendants. One can read from (14) the degree 
distribution of the root: 



1 



N(i-iy. 



N U ...,N,_ 1 



+ -+N l _ 1 ,N-l 



(22) 



Going over to the grand-canonical ensemble one finds: 

which, taking (16) into account and after integration 
yields 



(24) 



Using similar arguments one writes the weight of graphs 
where the root has the degree I and its daughter the 
degree k as 



z k i(N) 



Ql 



N(l - 2] 



^ATi+-+JV,_i,iV-lX 



JVi,...,JVl- 



xYlZ Ni z k {N t 



(25) 



i-i 



Hence 

Integrating the above equation one finally obtains 



Z k i(n(Z)) 



Ql 



qk 



(Z-2)!(fc-l)! 

I 



Z x l-2 r x 2 



o F{x 2 ) Jo 



dxi 



.k-l 



F{ Xl ) 



(27) 



which is the conditional probability, up to normalization, 
that a descendant has the degree k when the ancestor's 
degree is I. The normalization is determined summing 
over k on the r.h.s. above, with the result (/ — l)Zi(fi). 
It is a slightly different measure of neighbor correlation 
than that proposed in ref. [6], but it carries similar in- 
formation. Because of the integration the dependence on 
k and I does not factorize, in general. 

Eq. (27) holds in the grand-canonical ensemble. One 
would like to have an expression valid in the canonical 
ensemble, where the graph has a well defined number of 
nodes. This requires Laplace transforming the r.h.s. of 
eq. (27) (one should transform first and normalize next). 
One can argue that in the thermodynamical limit and 
in the regime 5fi — > this often amounts to just replace 
Z -> x in (27). Indeed, let Z - x ~ 8fi a . Then 



d6n e 5 ^ N Zki(n) 



Zki(fi) J 



SfiN—c 8fi '~ 



(28) 



Evaluating the integral by the saddle-point method one 
gets an exponential of a quantity scaling like ./V"^" -1 ). 
This gives a factor of unity in the limit N — > oo, provided 
< a < 1, a condition often met in applications. 



■5 



D. Fractal dimension 

Repeating over and over the iteration process leading 
to eq. (27) one gets 



Zk 1 k 2 ...k r (f*(Z)) = 



n 

J=2 



{kj - 2)! (h - 1)! J Q F(x r ) 



dx r —!- — - x 



■F(*l) 



Summing over node degrees fci, k 2 , ■ ■ ■ ,k r one obtains the 
weight of all graphs with a point separated by r links from 
the root, ie the two-point correlation function C(r, 
introduced in sect. 1.2 : 



C(r,n{Z)) = 



Jo r F{x r ) J 



F'(Xr-l) 

r • • • x 



[ X3 dx 2 ^± r& xi (30) 

Jo F(x 2 ) Jo 



For finite x, replacing the upper limit of integration over 
x\ by x and performing all integrations, one gets 



C(r,n{Z))<x 



(mF(Z)) 
(r-1)! 



r-1 



(31) 



Hence, the tail of C(r, /j.) falls at least as fast as a Pois- 
sonian. Consequently (r)^ grows at most like hiF(Z). 
Assuming that F(z) has at most a power singularity at 
z = x one concludes that 



and therefore 



{r)p < const In -!- 



( r ) n < const In N 



(32) 



(33) 



since Sfi scales like A -1 . The argument is rather heuris- 
tic, but suggestive (see also the examples in the next sec- 
tion). It appears that generically the causal trees have 
the small-world property du = oo, contrary to the max- 
imum entropy trees whose generic fractal dimension is 
finite [3, 39, 41]. This phenomenon is easy to under- 
stand intuitively : the causal structure suppresses long 
branches. This can be seen by noting that along a branch 
from the root to the leaf no label permutations are pos- 
sible, hence a tree with a few long branches admits much 
less causal labelings then a "short fat" one. 



E. Synchronic view of growing networks 

So far, our discussion was very general. Let us now es- 
tablish a bridge to the popular growing network models. 



Of course, here we consider only those models, where one 
constructs tree graphs. Successive nodes are attached, 
one at a time, the attachment probability being a func- 
tion of the degree of the target node : 



Prob(/c | t) 



A(t) 



ith A(t) = Y J N k it)A k (34) 



where Nk(t) is the number of vertices with degree k in 
the tree at time t. It is obvious that the tree constructed 
that way has a causal structure: nodes are labeled by the 
attachment time. 

Notice, that from the perspective of a model builder, 
the concepts of causal and growing networks are com- 
plementary rather than equivalent. A causal network 
is defined by specifying the microstate weights. Every 
growing network is causal, of course, but the weights 
corresponding to a given growth process can have a very 
complicated, non-local structure. And conversely, given 
a set of weights it is, in general, not evident what is the 
corresponding growth process. Only for a class of models 
there exists a stationary attachment kernel Ak- 

For linear or shifted linear attachment kernels A k the 
normalization factor 1/A(t) depends only on the size of 
the previous configuration and is therefore the same for 
all trees of the same size [37]. Hence, the preferential 
attachment recipe is compatible with the factorization of 
p(T): Working in the canonical ensemble, at fixed time, 
we can drop the normalization factor altogether without 
any loss of generality and set [38] 



In = Qi Yi A ^ n > 1 

fe=i 



(35) 



where q\ is some positive constant (eventually set to 1 
in explicit calculations). It is instructive to check by in- 
spection that the graph weights produced by the recur- 
sion relation (14), coincide with those generated by the 
growing network model recipe. 



F. Examples 

1. Barabasi- Albert model [4] 

In this model A k = k. Therefore q n = (n— 1)!, F(x) — 
(1 — x)^ 1 and x = 1. Thus G(x) = x — ^x 2 and the 
solution to the equation (18) is 



Z(n) = 1 - Vl - 2e-A* rj 1 - Jb^x (36) 
The degree distribution is found from (20) : 

4 



tt„ = 2 f dx (1 - xf x™- 1 
Jo 



n(n + l)(n + 2) 
which coincides with the solution given in refs. [6, 36]. 



(37) 
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The near-neighbor correlation is readily found from 
(27) : 

Z k i = {I -I) dx 2 (1 - x 2 ) x l 2 - 2 / dx! (1 - xi) x*- 1 
Jo Jo 

k{k + l){k + 1 - l){k + 1) ^ 

This is very similar to the result of ref. [6], where a 
slightly different quantity has been calculated. The phys- 
ical content is the same: the causal structure has induced 
correlations between the node degrees. 
Since in this example 



3. Constant attachment kernel 



s>ms>-s: 



FMJ a dn = /„ ''"If©- 1 ! < 39) 



one has 



,r-l 



C(r, (i(Z)) = {ln [ iZ) i -C(r-1, (i(Z)) (40) 



Hence 



(r-1)! 



(r>„~]nF(Z)~ln 



which implies the following scaling law : 



(41) 



(42) 



For A k = 1 



^) = EfW^ 1=eI 



(48) 



The interest of this example is in the infinite radius of 
convergence of the above series: x = oo. The degree 
distribution is 



OC 



dxe Zx x 



2x^71— 1 c\ — n 



(49) 



a result found by several people, including the authors of 
ref. [6] . We again skip the calculation of Z^i ■ 
Since d\nF(x) — dx one finds 



and using 



one derives 



C(r,(i{Z)) = 



Z = In — , Sfj, -» 

0(1 



(r) N ~lnN 



(50) 



(51) 



(52) 



2. Krapivsky-Redner model [6] 

Now Ak = k + w, q n = r(n + w) and 

F{x) =T{w + l){l-x)- w -\ 
Of course x = 1 and 



G(l) = 



1 



(2 + w)I> + 1) 



(43) 



(44) 



Using the above and evaluating the Euler integral that 
appears in the present case on the r.h.s. of (20) one 
obtains 



(2 + w)r(3 + 2w) T(n + w) 
T(l + w) r(3 + n + 2w) 



(45) 



Again this reproduces exactly the result of ref. [6] . The 
reader can easily calculate the correlation Z^i- We skip 
this calculation here, because the result is not particu- 
larly instructive: the remarkable fact is the very existence 
of the correlation, not its particular form, which in this 
particular case is rather cumbersome. 

The calculation of C(r,(i) is identical to that carried 
out for the Barabasi- Albert model, except that 



Z=l- 5(i— ,6(i^Q 



which implies 



(r) 



N 



w - 
w 



■In AT 



(46) 



(47) 



III. CAUSAL VERSUS MAXIMALLY RANDOM 
TREES 

A. An important tiny difference 

It is instructive to consider for also the well known case 
of maximum entropy trees (cf. refs. [3, 39-41]). This 
will help putting the results of the preceding section in 
proper perspective. As in the preceding section we start 
with the composition rule for trees, which in the present 
case, when the causality constraint is lifted, reads 



£(Ti 



•T k ) 



(N + iy. i 

N!\...N k \k\ 



L(T 1 )---L(T k ) (53) 



Compared to (11) the difference may seem tiny: one just 
has (N+ 1)! instead of AH in the numerator on the r.h.s.. 
It is so because there is no causality constraint and there- 
fore all N + 1 labels can be arbitrarily distributed among 
the trees. However, this tiny difference has rather dra- 
matic consequences. Indeed, repeating the steps which 
led from (11) to (16) one obtains 



Z{p) - e^F(Z) 



(54) 



No derivative appears on the left-hand side! Instead of 
eq. (18) one has 



e-» = H(Z) = 



F{Z) 



(55) 
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where the function H(Z) plays the role analogous to that 
of G(Z) in the preceding section. Eq. (55) implies that 
the vertex degree distribution is 



n-l 



1 gnZ* 

F(x t ) (n-l)! 



(56) 



Here x* = min(x max , x) and a; TOO x is the position of 
the maximum of the function H(x) : H'{x max ) — 
0, H"{x max ) < 0. 

Since no derivatives appear as one goes over to the 
grand-canonical ensemble, one does not integrate either. 
Therefore, iterating the recursion relation one finds that 
Zki (fj) and C (r, jj) factorize. The latter equals 



C(r,u)cx {e-»F'{Z)) r 



(57) 



as first derived by Ambj0rn et al [39] using diagrammatic 
arguments. 

Let us now consider three examples, with q n chosen as 
in sects. IIF.l - IFF. 3 respectively. 



B. Examples 

1. First example: q n = (n — 1)! 

With this choice, that of sect. IIF.l, the function H(x) 
reads H(x) = x — x 2 and has its maximum at x = 1/2. 
Hence x* = 1/2. Since F(x») = 2, the degree distribution 
calculated from (56) is now 

TT n = 2"" (58) 

as in sect. II.F.3 ! 

It is easy to check that for small 5(i one has 

e-»F'(Z(p)) = 1 - y/sjl (59) 

and therefore 

C(r, M )cxe-^ r (60) 

which implies 

(r) N ~ iV 1/2 (61) 
Hence, dn = 2, the generic value [39]. 

2. Second example: q n = F(n + w) 

With the choice of sect. IIF.2 the function H(x) = 
T(w + l)x(l — x) w+1 has its maximum at x max = 1/(2 + 
w) , which is smaller than the radius of convergence x — 1 
of the series F(x) as long asw > —1. Thus x* = l/(2+w) 
and 



2 + w /l + w\ 1+w T(n + w) 
T(l + w) \2 + w) T(n) 



1 



2 + w 



(62) 

The fall of the degree distribution is again exponential. 
One easily checks that again dn = 2. 



3. Third example: q n = 1 

With the choice of sect. IIF.3 the function H(x) = 
xe~ x has its maximum at x max = 1- Thus x* = 1 and 
correspondingly 



(63) 



(n-l)! 



One can easily check that this result holds for a more 
general family of weights q n = a/3™ -1 independently of 
the values of a and /3, as long as they are positive. Again 
d H = 2. 



C. How to get identical degree distribution in the 
two models 

It is easy to adjust the input parameters of the two 
models to obtain identical degree distributions: Suppose 
that in the causal model the degree distribution is n n . 
For a connected tree one necessarily has 



^2 mr n = 2 

n 

Set in the maximally random model 



Obviously 



q n = (n - 1)! tt„ 



[ff- 1 ]'(l)^^(n-2) 7 r„ = 



(64) 



(65) 



(66) 



and [77 _1 ]"(1) > 0. Thus, H(x) has a maximum at 
x = 1 . One easily convinces oneself that this is the only 
maximum of that function. Furthermore, by assump- 
tion (see (64)), the radius of convergence in F(x) is > 1. 
Hence x* = 1 and 



random 



(67) 



by virtue of eq. (56). However, the coincidence of the 
degree distributions does not imply that geometries are 
similar. On the contrary, as already emphasized, the 
graph linear sizes have generically a very different scal- 
ing behavior in the two models (see also some computer 
simulation results in [42]). 



IV. SUMMARY AND CONCLUSION 

The results of this paper are another illustration of 
the claim that the opposition between diachrony and syn- 
chrony is to large extent an illusion, except if one is inter- 
ested in very specific phenomena, like aging, intrinsically 
reflecting the running of time. We have discussed the 
geometry of networks endowed with a causal structure 
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using the conventional framework of equilibrium statisti- 
cal mechanics. Hence, models that are usually described 
by specific master equations and static, maximum en- 
tropy models have been treated alike. We focused on tree 
graphs, because only for trees we are able to proceed an- 
alytically. We have derived general formulae describing 
the degree distribution, the ancestor-descendant correla- 
tion and the probability a node lives at a given geodesic 
distance from the root. Using this last results have shown 
that our causal networks have generically the small- world 
property, it their Hausdorff dimension is infinite. 

We have also compared the causal model with the 



maximally random one, assuming the same microstate 
weights. Because of different symmetry properties of the 
graphs - in the causal model only a subclass of labelings is 
allowed - the geometries are dramatically different: in the 
causal model the degree distribution is qualitatively dif- 
ferent, inter-node correlations are induced and the Haus- 
dorff dimension becomes infinite instead of being finite. 
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